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Investigation  of  the  role  in  which  the  microporous  layer  (MPL)  affects  the  overall  diffusion  of  gases  and 
conduction  of  heat  in  polymer  electrolyte  membrane  (PEM)  fuel  cells  is  of  great  interest.  In  this  study, 
we  used  stochastic  models  to  generate  a  three-dimensional  reconstruction  of  the  MPL.  This  work  is  a 
continuation  of  the  methods  presented  in  [10].  A  parametric  study  was  carried  out  to  investigate  the 
effects  of  the  MPL  structure  and  MPL  porosity  on  its  diffusion  coefficient  and  thermal  conductivity.  It  was 
found  that  increasing  the  volume  of  the  small  pores  of  the  MPL  while  keeping  its  overall  porosity  constant 
results  in  an  increase  in  the  Knudsen  diffusion;  hence  a  decrease  of  the  overall  diffusion  coefficient.  This 
similar  trend  is  observed  again  once  the  porosity  of  the  MPL  is  decreased.  An  increase  in  the  volume  of 
small  pores  also  resulted  in  an  increase  of  the  thermal  conductivity  of  the  MPL.  The  parametric  study  was 
also  extended  to  understand  the  effect  of  applying  the  MPL  onto  the  gas  diffusion  layer  (GDL).  In  this  case, 
we  investigated  the  effect  of  MPL  thickness,  porosity  and  its  penetration  into  the  GDL.  The  effect  of  the 
thickness  on  the  thermal  conductivity  and  diffusion  coefficient  of  an  MPL/GDL  assembly  can  be  explained 
using  a  resistance  network.  An  increase  in  the  penetration  depth  of  the  MPL  results  in  an  increase  of  the 
thermal  conductivity  and  a  decrease  of  the  diffusion  coefficient. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  microporous  layer  (MPL)  is  a  layer  composed  of  carbon 
agglomerates,  which  is  applied  to  the  interface  between  the  gas 
diffusion  layer  (GDL)  and  the  catalyst  layer  (CL)  of  polymer  elec¬ 
trolyte  membrane  (PEM)  fuel  cells.  This  layer  is  used  to  achieve  a 
smooth,  continuous  interface  between  the  GDL  and  the  CL.  The  MPL 
consists  of  a  powdery  mixture  of  carbon  black  and  PTFE  particles, 
is  sprayed  onto  the  GDL  and  then  sintered  so  the  PTFE  can  bind 
the  powder  together.  Experimental  and  numerical  investigations 
have  revealed  that  the  performance  of  the  cell  is  improved  with 
the  use  of  the  MPL  and  much  work  is  dedicated  to  understand  the 
mechanisms  in  which  the  MPL  enhances  the  performance. 

One  of  the  main  mechanisms  in  which  the  MPL  affects  the  over¬ 
all  performance  of  the  cell  is  through  assisting  with  liquid  water 
management.  This  phenomenon  has  been  studied  numerically  in 
[1,2],  where  in  both  studies,  a  two  phase  model  was  used  to  inves¬ 
tigate  the  effects  of  placing  a  microporous  layer  at  the  cathode  side. 
The  results  indicate  that  the  MPL  reduces  the  overall  water  flux 
to  the  cathode  GDL  particularly  when  the  liquid  water  saturations 
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are  high;  hence,  improving  membrane  hydration  and  reducing  the 
probability  of  GDL  flooding.  The  reduction  in  the  water  flux  towards 
the  cathode  side  is  mainly  due  to  the  increased  tendency  of  the 
water  to  flow  towards  the  anode.  In  other  words,  the  presence  of  an 
MPL  at  the  cathode  side  helps  in  establishing  another  route  for  liq¬ 
uid  water  removal  from  the  cell.  However,  this  trend  is  not  observed 
when  placing  the  MPL  at  the  anode  side  only  as  Weber  and  Newman 
reported  in  [1].  From  their  findings,  they  suggested  that  the  MPL 
improves  the  performance  of  the  cell  due  to  the  following  reasons: 


1.  The  addition  of  the  microporous  layer  reduces  the  ohmic  losses 
in  the  fuel  cell.  This  is  mainly  attributed  to  the  enhanced  hydra¬ 
tion  of  the  electrolyte  membrane  and  the  reduction  of  the  overall 
resistance  between  the  GDL  and  the  CL. 

2.  The  MPL  reduces  liquid  water  flooding  in  the  GDL;  hence,  the 
transport  of  gaseous  species  is  enhanced  and  the  mass-transport 
limitations  are  minimized. 

3.  The  MPL  is  a  region  that  enhances  the  chemical  and  mechanical 
stability  of  the  CL  and  the  membrane.  During  compression,  with¬ 
out  the  MPL,  the  catalyst  layer  may  become  entrenched  into  the 
GDL;  hence,  hindering  its  overall  activity  and  accessibility.  Sim¬ 
ilarly,  during  compression,  the  GDL  may  pierce  the  membrane 
resulting  in  reactant  crossover. 
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Fig.  1.  Pore  distribution  of  various  microporous  layers  -  data  extracted  from  [3-6]. 


1.1.  Structure 

The  structure  of  the  MPL  is  complicated  in  that  the  distribution 
of  the  carbon  agglomerates  highly  depends  on  the  manufacturing 
and  heat  treatment  processes,  which  directly  affect  the  pore  dis¬ 
tribution  of  the  MPL  as  is  evident  by  Fig.  1.  The  data  presented  in 
this  figure  has  been  extracted  from  [3-6].  Kumbur  et  al.  [3,4]  eval¬ 
uated  the  pore  distribution  for  SGL  24DC  and  SGL  24BC.  One  of  the 
main  contributors  to  the  differences  between  these  two  SGL  car¬ 
bon  papers  could  be  associated  with  the  amount  of  PTFE  used  to 
treat  the  MPL.  SGL  24DC  contains  more  PTFE  than  SGL  24BC.  The 
MPL  measured  by  Martinez-Rodriguez  et  al.  [5]  was  measured  for 
an  MPL  applied  to  a  TORAY  TGP-H-060  with  10  wt.%  PTFE.  Finally, 
Williams  et  al.  [6]  investigated  the  pore  distribution  of  SGL  10BB 
carbon  paper.  Although  the  MPL  measured  in  [6]  and  in  [3,4]  are 
both  SGL  carbon  papers,  their  manufacturing  process  may  differ 
and  their  total  thicknesses  differ.  The  differences  in  weight  and 
thickness  of  the  carbon  paper  are  of  significance  since  the  pore 
volume  presented  in  all  the  studies  is  normalized  by  these  two 
properties. 

The  effect  of  the  manufacturing  and  application  processes  on  the 
structure  of  the  MPL  onto  the  carbon  paper  GDL  has  been  investi¬ 
gated  by  Chen  et  al.  [7].  In  their  study,  they  examined  the  effect 
of  two  preparation  methods,  dry  and  wet  layering  of  the  MPL.  The 
main  differences  between  these  two  processes  are  the  preparation 
of  the  carbon  composite  paste  and  the  heat  treatment  process.  From 
the  findings  of  their  investigation,  they  reported  visual  differences 
in  the  resultant  MPL,  which  in  turn  significantly  affected  the  overall 
performance  of  the  cell. 

It  is  clear  that  the  impact  of  the  MPL  on  the  overall  performance 
of  the  cell  is  positive.  Flowever,  more  insight  into  the  mechanisms 
behind  this  effect  is  still  needed  and  numerical  simulation  is  often 
used  to  visualize  the  MPL’s  local  impact.  To  do  so,  accurate  trans¬ 
port  properties  are  crucial  for  the  numerical  simulation  in  order  to 
ensure  its  accuracy.  Lienee,  estimation  of  the  transport  properties 
of  the  MPL  have  been  the  focus  of  many  studies  as  discussed  in  the 
following  section. 

1.2.  Properties  of  the  MPL 
1.2.1.  Diffusion  coefficient 

The  diffusion  coefficient  of  the  microporous  layer  is  essential 
for  simulating  gas  transport  in  the  layer  and  experimental  proper¬ 
ties  of  this  property  can  be  found  in  literature.  As  the  microporous 


layer  cannot  be  found  as  a  stand-alone  layer,  it  is  applied  to  one 
side  of  the  carbon  paper  GDL  and  measurements  are  made  for  an 
MPL/GDL  assembly.  Thus,  the  change  in  diffusion  coefficient  due 
to  the  addition  of  an  MPL  is  often  reported  in  literature.  LaManna 
and  Kandlikar  [8],  measured  the  effect  of  the  MPL  on  three  SGL 
carbon  paper  electrodes.  They  found  that  the  presence  of  the  MPL 
reduces  the  overall  through-plane  diffusion  coefficient  of  the  elec¬ 
trode  by  39%.  Chan  et  al.  [9]  also  measured  the  effect  of  the  MPL  on 
the  through-plane  diffusion  coefficient  of  SolviCore  carbon  paper 
using  a  Loschmidt  cell  apparatus.  In  their  study,  they  used  a  resis¬ 
tance  network  analysis  to  obtain  the  diffusion  coefficient  of  an 
oxygen-nitrogen  mixture  through  the  microporous  layer  and  esti¬ 
mated  it  as  1 .54  ±  0.22  x  1 0-6  m2  s-1  for  a  calculated  porosity  value 
of  64%. 

Numerically,  Becker  et  al.  [10]  developed  a  mathematical  for¬ 
mulation  to  take  into  account  the  Knudsen  diffusion  in  the  small 
pores  of  the  MPL.  They  showed  that  in  the  presence  of  the  MPL, 
the  diffusion  coefficient  in  both  the  in-  and  through-plane  direc¬ 
tions  decreases.  The  decrease  in  the  in-plane  direction,  however, 
was  not  as  pronounced.  Further,  they  found  that  the  change  in  the 
diffusion  coefficient  is  highly  dependent  on  the  overall  penetration 
of  the  MPL  into  the  carbon  paper  GDL.  Numerical  estimations  of 
this  nature  are  of  importance  since  much  about  the  structure  of  the 
MPL  is  unknown  or  is  very  difficult  to  measure. 

1.2.2.  Thermal  conductivity 

Estimation  of  the  thermal  conductivity  of  the  microporous  layer 
in  the  literature  has  recently  gained  interest.  Using  the  guarded 
heat  method,  Karimi  et  al.  [11]  reported  the  through-plane  ther¬ 
mal  conductivity  of  the  MPL  at  70  °C  and  for  various  compression 
pressures.  For  SolviCore  carbon  paper  with  an  MPL  and  30%  wet 
proofing,  they  found  the  thermal  conductivity  to  lay  in  the  range 
0.25-0.52  WK-1  m  for  a  compression  force  of  0.7-13.8 bar.  The 
in-plane  thermal  conductivity  of  this  SolviCore  paper  was  later 
measured  by  the  same  group  [12]  and  reported  as  3.87  WK-1  m. 
Burheim  et  al.  [13]  measured  the  through-plane  thermal  con¬ 
ductivity  of  an  MPL-electrolyte  membrane  assembly.  From  these 
measurements,  they  estimated  the  thermal  conductivity  of  the  MPL 
to  be  0.6  WI<-1  m.  The  thermal  conductivity  of  the  MPL  was  also 
measured  by  Unsworth  et  al.  [14].  They  used  the  measured  ther¬ 
mal  conductivity  of  a  bare  carbon  paper  and  a  carbon  paper/MPL 
assembly  to  isolate  the  contribution  of  the  MPL.  They  estimated 
the  thermal  conductivity  of  the  MPL  to  be  0.30  WIG1  m.  They  also 
found  that  the  effect  of  the  MPL  on  the  through-plane  thermal 
conductivity  of  the  electrode  highly  depends  on  the  amount  of 
compression  pressure  applied  to  the  electrode.  As  the  compression 
pressure  increases,  the  effect  of  the  MPL  on  the  thermal  conductiv¬ 
ity  becomes  negative.  With  this  finding,  they  concluded  that  unlike 
carbon  paper,  the  structure  of  the  MPL  does  not  change  due  to  com¬ 
pression.  Hence,  the  overall  thermal  conductivity  of  the  MPL  does 
not  change  with  compression.  Under  compression,  however,  the 
contact  between  the  carbon  fibers  in  the  carbon  paper  increases 
and  its  thermal  conductivity  increases;  thus,  at  these  conditions, 
the  MPL  acts  as  a  thermal  resistor  and  decreases  the  overall  thermal 
conductivity  of  the  backing  layer. 

In  their  study,  Burlatsky  et  al.  [15]  reported  their  internal  mea¬ 
surements  of  the  thermal  conductivity  of  the  MPL  as  a  stand  alone 
layer.  The  measurements  reported  show  that  the  thermal  conduc¬ 
tivity  of  an  MPL/GDL  assembly  is  much  lower  than  that  of  a  GDL 
only.  They  also  showed  that  the  measurements  varied  depending 
on  the  measurement  technique.  For  instance,  the  thermal  conduc¬ 
tivity  of  a  thin  MPL  with  low  and  high  Teflon  content  was  0.035 
and  0.057  WK_1m,  respectively,  using  the  laser  flash  method. 
However,  with  the  modified  substrate  method,  the  value  reported 
was  not  affected  by  the  Teflon  content  and  it  was  reported  to  be 
0.097  WK-1  m.  Comparing  the  two  methods,  one  can  see  there 
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is  much  variation  depending  on  the  method.  Specific  conclusions 
about  this  study  cannot  be  drawn  as  much  information  surrounding 
the  experimental  apparatus  and  experimental  conditions  is  uncer¬ 
tain. 

The  thermal  conductivity  of  the  MPL  was  also  reported  theoret¬ 
ically  in  literature.  In  their  model,  Kang  and  Ju  [16]  analyzed  the 
effects  of  the  microporous  layer  on  the  transport  mechanisms  in 
PEM  fuel  cells.  In  their  study,  they  assumed  the  thermal  conductiv¬ 
ity  of  the  MPL  to  be  1 .0  W  K-1  m. 

As  can  be  seen,  published  data  of  the  thermal  conductivity  of  the 
MPL  vary  tremendously  and  are  affected  by  the  method  of  measure¬ 
ment.  This  variation  coupled  with  limited  knowledge  in  regards  to 
the  carbon  used  in  the  manufacture  of  the  MPL  makes  it  difficult  to 
judge  its  thermal  conductivity.  Hence,  in  this  study,  we  investigate 
the  effect  of  the  agglomerate  structure  on  the  thermal  conductivity. 

1.2.3.  Electrical  conductivity 

Similarly  to  the  diffusion  coefficient  and  the  thermal  conduc¬ 
tivity,  the  electrical  conductivity  of  the  MPL  is  measured  based  on 
the  change  in  the  overall  electrical  conductivity  of  the  MPL/carbon 
paper  assembly  in  comparison  to  that  of  the  bare  carbon  paper. 
Experimental  measurements  of  the  electrical  conductivity  of  the 
backing  layer  in  the  presence  of  the  MPL  are  scarce.  Ismail  et  al. 
[17]  reported  the  contact  resistance  between  the  porous  layer  and 
the  test  rig  in  the  presence  of  the  MPL.  They  showed  that  the  con¬ 
tact  electrical  resistance  in  the  presence  of  the  MPL  decreases  with 
the  highest  change  observed  was  an  increase  by  40%  in  the  overall 
contact  resistance.  The  change  in  the  in-plane  electrical  conductiv¬ 
ity  in  the  presence  of  the  MPL  was  reported  by  Williams  et  al.  [6] 
where  they  found  the  in-plane  conductivity  to  increase  by  4%.  This 
small  change  in  comparison  to  its  through-plane  counterpart  sug¬ 
gests  that  the  conduction  of  the  electrons  in  the  in-plane  direction 
via  the  carbon  fibers  is  efficient. 

1.2.4.  Intrinsic  permeability 

The  intrinsic  and  relative  permeability  can  be  used  as  another 
measure  to  understand  the  effect  of  the  MPL  on  the  transport  of 
the  gas  in  the  electrode.  Again,  due  to  the  difficulty  of  separating 
the  microporous  layer  from  the  carbon  paper,  the  permeability  of 
the  bare  carbon  paper  is  first  measured  and  then  the  permeabil¬ 
ity  of  the  carbon  paper/MPL  assembly  is  measured.  Hence,  one 
can  obtain  an  estimation  of  the  effect  of  the  MPL  on  this  prop¬ 
erty.  The  through-plane  intrinsic  permeability  of  SGL  10BA  (bare) 
and  SGL  10BB  (with  MPL)  was  measured  by  Williams  et  al.  [6].  In 
their  study,  they  reported  the  intrinsic  permeability  to  decrease 
from  3.1  x  1 0-1 1  m2  to  5.8  x  1 0-13  m2 ;  a  decrease  of  about  2  orders 
of  magnitude.  The  hinderance  of  the  permeability  in  the  through- 
plane  direction  is  mainly  attributed  to  the  small  pores  present  in  the 
MPL  and  the  introduction  of  the  Knudsen  flow.  The  same  trend  is 
observed  by  Ihonen  et  al.  [  1 8  ]  where  they  reported  the  through-  and 
in-plane  intrinsic  permeability  of  SGL  10BA  (bare)  and  SGL  10BC 
(with  MPL).  The  in-  and  through-plane  permeability  of  SGL  10BA 
was  measured  as  3.3  x  10-11  and  1.8  x  10-11  m2,  respectively  and 
the  in-  and  through-plane  permeability  of  SGL  10BC  was  measured 
as  2.2  x  10-11  and  3.3  x  10-13  m2,  respectively.  As  expected,  the 
change  in  the  in-plane  permeability  due  to  the  presence  of  the  MPL 
is  not  as  pronounced  as  its  through-plane  counterpart,  which  is 
mainly  due  to  the  availability  of  the  pore  space  in  the  carbon  paper 
GDL  in  the  in-plane  direction  for  fluid  transport. 

1.2.5.  Capillary  pressure 

The  capillary  pressure  is  a  measure  for  understanding  liquid 
water  transport  in  the  porous  material.  Using  the  method  of  stan¬ 
dard  porosimetry,  Kumbur  et  al.  [3,4,19]  measured  the  capillary 
pressure  for  SGL  24  series  with  an  MPL  and  proposed  a  correla¬ 
tion  for  the  capillary  pressure  based  on  the  compression  level,  PTFE 


content  and  temperature.  However,  in  their  work,  they  did  not  sep¬ 
arate  the  effect  of  the  MPL  from  that  of  the  carbon  paper  GDL  as 
they  did  not  measure  the  capillary  pressure  for  a  bare  carbon  paper. 
Hence,  a  full  understanding  of  the  effect  of  the  MPL  on  the  overall 
capillary  pressure  could  not  be  obtained.  In  the  study  by  Gostick 
et  al.  [20],  the  authors  measured  the  capillary  pressure  of  SGL  1 OBA 
(bare)  and  SGL  10BB  (with  MPL)  and  were  able  to  isolate  the  effect 
of  the  MPL  on  the  capillary  pressure.  Their  finding  is  interesting  as 
it  suggests  that  the  MPL  plays  a  significant  role  in  controlling  liquid 
water  levels  in  the  GDL.  They  used  this  finding  to  further  analyze 
the  effect  of  the  MPL  on  the  capillary  pressure  in  [21  ].  They  showed 
that  the  liquid  water  saturation  at  the  breakthrough  pressure  drops 
by  3%  in  the  presence  of  the  MPL.  This  suggests  that  a  lower  liquid 
water  saturation  level  would  be  present  in  the  porous  electrode 
and  the  MPL  plays  a  crucial  role  in  water  management  of  the  cell. 

As  can  be  seen,  much  experimental  work  has  been  done  to  esti¬ 
mate  the  properties  of  the  MPL.  However,  due  to  the  uncertainty 
related  to  the  structure  of  the  MPL,  it  is  difficult  to  obtain  a  com¬ 
prehensive  understanding  of  these  properties.  In  this  study,  we 
estimate  the  diffusion  coefficient  and  thermal  conductivity  of  the 
microporous  layer  numerically  using  a  three-dimensional  recon¬ 
struction  of  the  layer.  We  examine  the  effect  of  the  following 
parameters  on  the  transport  properties: 


•  Effect  of  MPL  porosity  -  experimentally,  it  is  difficult  to  accurately 
measure  the  porosity  of  the  MPL  due  to  its  penetration  into  the 
GDL.  Efforts  have  been  made  to  estimate  this  value  as  seen  in 
[9,20]. 

•  Effect  of  MPL  thickness  -  in  reality,  the  thickness  of  the  MPL 
varies  depending  on  the  sample  under  study  as  has  recently  been 
reported  by  [9]. 

•  Effect  of  MPL  penetration  -  the  penetration  thickness  of  the  MPL 
has  a  direct  effect  on  the  overall  porosity  of  the  electrode;  hence, 
it  affects  the  contact  area  between  the  carbon  fibers  of  the  GDL. 


2.  Structures 

For  common  GDL  and  MPL  materials,  the  pore  sizes  of  MPL  and 
GDL  differ  by  more  than  two  orders  of  magnitude.  This  implies,  that 
a  pore-scale  model  (with  a  number  of  grid  points  manageable  by  the 
subsequent  numerical  simulations)  cannot  at  the  same  time  show 
a  representative  part  of  the  GDL  and  resolve  the  pores  of  the  MPL. 
Therefore,  we  follow  the  approach  taken  in  [10]  and  create  two 
different  structure  models.  First,  a  high-resolution  model  which 
only  shows  the  pore  structure  of  the  MPL.  Second,  a  low-resolution 
model  showing  the  whole  thickness  of  the  GDL,  where  the  MPL  is 
added  as  homogenized  porous  medium. 

In  both  cases,  the  structure  models  are  created  by  stochastic 
methods.  This  means,  that  fibers  or  particles  are  distributed  ran¬ 
domly  in  a  predefined  volume.  These  objects  are  then  used  to  create 
the  voxel  mesh,  which  serves  as  mesh  for  the  numerical  simula¬ 
tions. 


2.1.  MPL 

The  three  dimensional  structure  of  the  reconstructed  MPL  is 
given  in  Fig.  2.  Scanning  electron  microscope  (SEM)  images  (see  e.g. 
Fig.  3b  or  [10,22])  show  that  the  MPL  consists  of  carbon  particles 
which  are  -  roughly  -  of  spherical  shape  and  equal  size.  These  par¬ 
ticles  are  not  evenly  distributed,  but  form  a  complicated  network 
of  large  pores  and  agglomerates.  This  observations  motivated  the 
following  construction  algorithm,  which  is  also  illustrated  in  Fig.  4: 
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•  In  a  first  step,  part  of  the  volume  is  set  aside  for  the  large  pores. 
This  is  done  by  randomly  distributing  large  spheres,  which  are 
allowed  to  overlap. 

•  In  the  second  step,  the  carbon  particles  are  created  by  distribut¬ 
ing  small  spheres,  which  are  again  allowed  to  overlap  with  each 
other.  In  this  step,  the  centers  of  the  small  spheres  are  not  allowed 
to  be  inside  one  of  the  large  spheres  of  the  previous  step. 

•  In  the  third  and  last  step,  the  small  particles  are  glued  together 
by  filling  the  small  pores  in  between  them. 

This  basically  allows  to  change  5  parameters:  In  the  first  step,  the 
diameter  dpore  of  the  spheres  and  the  volume  fraction  r]pove  filled  by 
them.  In  the  second  step,  the  diameter  dparticle  °f  the  spheres  and 
the  volume  fraction  particles  filled  by  them.  In  the  third  step,  the 
maximum  pore  diameter  d^\\  to  be  filled. 

In  this  study,  the  parameters  dparticle  =  40  nm  and  dfill  =  20nm 
were  kept  fixed  while  the  others  were  varied.  All  MPL  models  used 
had  a  domain  size  of  (3  pan)3,  consisting  of  6003  cubic  voxels  with 
a  voxel  length  of  5  nm. 


2.2.  GDI 


Fig.  3.  (a)  Cross-section  of  the  reconstructed  microporous  layer;  (b)  SEM  image  of 
a  microporous  layer  [22]. 


Following  the  ideas  in  [10,23],  a  GDL  similar  to  a  Toray  paper 
can  be  modeled  by  straight,  unbounded  fibers  with  a  circular  cross 
section  of  7  pan.  For  the  GDL  model  used  in  this  study,  the  fiber 
positions  and  directions  were  found  by  a  Poisson  line  process  with 
an  orientation  tensor  of 

/  0.499975  0  0  \ 

0  0.499975  0  (1) 

V  0  0  0.00005  J 

i.e.  fibers  are  mainly  oriented  in  the  xy-plane  (in-plane).  Fibers  were 
added  to  the  structure  until  a  porosity  of  90%  was  reached.  As  gas 
diffusion  layers  contain  both  fibers  and  binder,  50  wt.%  binder  were 
added  to  the  structure  using  the  method  described  in  [24],  where 
it  is  assumed  that  the  binder  basically  behaves  like  a  wetting  fluid 
and  fills  the  smaller  pores.  The  resulting  GDL  models  thus  have  an 
overall  porosity  of  80%. 

The  GDL  models  used  in  this  study  had  a  domain  size  of 
512  pm  x  512  pm  x  256  p,m,  consisting  of  512  x  512  x  256  cubic 
voxels  with  a  voxel  length  of  1 .0  pm. 


2.3.  GDL/MPL  assembly 

The  MPL  is  added  to  the  GDL  model  as  a  homogenized  porous 
layer.  Thus,  a  voxel  in  the  resulting  mesh  represents  either  solid 
(fibers  or  binder),  void  (GDL  pores)  or  porous  material  (MPL)  and 
the  computational  domain  Q  subdivides  into 


if2  —  ^solid  U  if2void  U  if2porous*  (2) 


Fig.  5  illustrates  one  of  the  resulting  structure  models.  The  com¬ 
bined  GDL/MPL  models  have  the  same  resolution/voxel  length  as 
the  GDL  models.  The  domain  size  is  512  pm  x  512  pm  in  the  xy- 
plane  and  256  p,m  plus  MPL  thickness  in  the  through-plane  (z) 
direction.  In  this  study,  we  vary  the  MPL  thickness  and  penetration 
depth. 
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Fig.  4.  Illustration  of  the  MPL  reconstruction  algorithm. 


3.  Property  calculation 

The  generated  3D  voxel  meshes  are  used  to  calculate  the  effec¬ 
tive  properties  of  the  porous  layer.  Here,  we  consider  pore  size 
distribution,  thermal  conductivity  and  diffusivity. 

At  first,  we  determine  the  properties  of  the  MPL  using  an  MPL 
structure  model.  The  resulting  properties  then  enter  into  the  sim¬ 
ulation  of  the  GDL/MPL  assembly:  As  the  MPL  is  not  resolved  on 
this  scale,  each  voxel  in  X2porous  represents  a  porous  medium  with 
the  calculated  effective  properties.  Of  course,  different  MPL  models 


Fig.  5.  A  stack  of  the  MPL  and  GDL  (a)  cross-section  view;  (b)  3-D  view  -  MPL  (yel¬ 
low),  binder  (green),  fibers  (red),  pores  (white).  (For  interpretation  of  the  references 
to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  the  article.) 

result  in  different  MPL  properties  which  then  will  lead  to  different 
properties  of  the  assembly. 

3.1.  Pore  size  distribution 

The  radius  of  a  pore  is  defined  by  embedding  spheres  into  the 
pore  volume :  a  voxel  belongs  to  a  pore  of  radius  r  if  r  is  the  radius  of 
the  largest  sphere,  embedded  in  £2pore,  containing  this  voxel.  In  this 
way,  a  radius  is  assigned  to  each  voxel  and  the  pore  size  distribution 
shows  the  distribution  of  these  radii. 

3.2.  Thermal  conductivity 

Heat  transport  through  a  porous  medium  is  macroscopically 
described  by  Fourier’s  law  of  conduction 

J  =  — feeffVT,  (3) 

where  j  is  the  heat  flux  and  VT  is  the  temperature  drop  over  the 
medium.  In  general,  the  medium  may  be  anisotropic  and  the  effec¬ 
tive  conductivity  /<eff  is  a  3  x  3  matrix.  Here,  the  MPL  models  are 
isotropic  by  construction  and  in  this  case,  Eq.  (3)  holds  for  a  scalar 

/<eff. 

The  conductivity  /<eff  can  be  determined  numerically  by  solving 
the  stationary  heat  equation 

V.(k(x)VT)  =  0  in£2  (4) 

directly  on  the  3D  structure  model.  Here,  /<(x)  is  the  local  (isotropic) 
thermal  conductivity,  which  was  set  to  k(x)  =  0  in  T2vo[d  and 
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/<(*)  = /<c arbon  in  ^soiid-  To  solve  (4),  the  EJ-Heat  solver  [26]  imple¬ 
mented  in  GeoDict  [25]  was  used. 

The  same  methodology  is  used  to  find  the  thermal  conductivity 
of  the  MPL/GDL  assembly.  In  this  case,  the  conduction  of  heat  occurs 
in  four  distinct  parts  of  the  assembly:  the  binder  of  the  GDL,  the 
solid  fraction  of  the  GDL,  the  void  region  of  the  GDL  and  the  MPL. 
Due  to  the  differences  in  the  scale  of  pores  of  the  two  layers,  the 
true  pore  morphology  of  the  MPL  cannot  be  coupled  with  that  of 
the  GDL  and  the  computed  thermal  conductivity  of  the  MPL  is  used. 
Mathematically  the  heat  conduction  is  represented  by  Eq.  (4)  with 
the  thermal  conductivity  k  varying  as: 

{0  X  G  ^void 

fempi  X  e  ^mpl  (5) 

^solid  X  G  *f^solid 

where  we  assume  that  the  binder  and  solid  fibers  of  the  GDL  are 
manufactured  with  the  same  material;  hence,  their  thermal  con¬ 
ductivity  is  equal. 

3.3.  Diffusion  coefficient 

In  general,  diffusion  occurs  due  to  collision  of  gas  molecules 
with  each  other  (bulk  diffusion)  and  to  collision  with  the  pore  walls 
(Knudsen  diffusion).  The  contribution  of  each  type  of  diffusion  to 
the  overall  transport  process  is  described  by  the  Knudsen  number 
Kn ,  which  is  defined  through 

Kn=f  (6) 


Fig.  6.  Effect  of  pore  density  and  diameter  on  the  overall  pore  distribution  of  the 
microporous  layer. 


we  determine  Dbulk.  To  solve  Eq.  (9)  numerically,  we  use  GeoDict 
[25]  and  the  underlying  explicit  jump  finite  volume  solver  [26].  The 
solver  works  directly  on  the  voxel  grid,  so  no  re-meshing  is  needed. 
This  enables  us  to  find  the  solution  fast  and  memory-efficient. 

To  determine  the  diffusivity  of  the  GDL/MPL  assembly,  we  have 
to  solve 


where  X  is  the  mean  free  path  and  lc  is  the  characteristic  length  of 
the  porous  media,  e.g.  the  average  pore  diameter.  Large  Knudsen 
numbers  Kn  » 1  indicate  that  the  pore  sizes  of  the  porous  medium 
are  much  smaller  than  the  mean  free  path  of  the  gas  and  diffusion 
occurs  in  the  form  of  Knudsen  diffusion.  Small  Knudsen  numbers 
Kn  <$c  1  indicate  that  the  pores  of  the  porous  medium  are  much 
larger  than  the  mean  free  path  of  the  gas  and  bulk  diffusion  occurs. 
In  the  intermediate  range,  both  mechanisms  have  to  be  taken  into 
account. 

As  the  mean  free  path  of  gases  is  in  the  order  of  1 00  nm,  diffusion 
inside  the  GDL  is  dominated  by  bulk  diffusion  and  the  Knudsen 
diffusion  can  be  neglected  in  this  case.  Inside  the  MPL,  pore  sizes 
are  of  the  same  order  of  magnitude  as  the  mean  free  path.  Therefore, 
both  diffusion  types  have  to  be  taken  into  account. 

To  determine  the  MPL  diffusivity,  we  follow  the  method  devel¬ 
oped  in  [10],  calculate  Knudsen  diffusivity  DI<n  and  bulk  diffusivity 
Dbuik  and  then  make  use  of  Bosanquet’s  formula  [27] 


-div(D(x)Vc)  =  0  in  £2  (11) 

instead  of  Eq.  (9),  where  the  local  diffusion  coefficient  D(x)  is  given 
as 

{t^void  X  G  ^void 

f^mpl  X  G  if2porous  (12) 

0  X  G  £2 solid 

Here,  Dmpl  is  the  previously  calculated  diffusivity  of  the  MPL  and 
DVOid  =  (1  /3)Xv  is  the  diffusion  coefficient  of  the  gas. 

As  mean  free  path  X  we  take  63.3  nm  and  as  mean  free  veloc¬ 
ity  v  =  425  m  s-1  [29].  To  solve  Eq.  (1 1 )  numerically,  we  again  use 
GeoDict  [25]  and  the  underlying  solver  [26]. 

4.  Results  and  discussion 

4. 1 .  Effect  of  MPL  structure 


D  =  (D 


-l 

bulk 


+  D 


-1)-1 
Kn  ' 


(7) 


to  determine  the  overall  diffusivity  D. 

We  determine  the  Knudsen  diffusivity  from  the  displacement  of 
the  diffusing  molecules: 


e  t 

D,(n  =  ^E[(xt  -X0 )(Xt  -Xo)  ] 


(8) 


Here,  xt-x0  is  the  displacement  of  a  molecule  at  time  t  from  its 
starting  position  x0, 6  denotes  the  porosity  and  E[.]  the  expectation 
value.  To  determine  the  displacement  of  each  molecule  we  apply 
the  random  walk  method  of  Babovsky  [28]  (also  described  in  [10]). 

We  find  the  bulk  diffusivity  with  the  classical  continuum 
mechanics  approach  by  solving  Laplace’s  equation 


-Ac  =  0  in  £2VO\d 


(9) 


In  this  section,  we  investigate  the  effect  of  the  MPL  structure  on 
its  diffusion  coefficient  and  thermal  conductivity.  Such  an  investi¬ 
gation  is  interesting  due  to  the  complications  associated  with  the 
reproducibility  of  the  structure  of  the  MPL.  In  order  to  carry  out  this 
comparison,  we  vary  the  pore  distribution  of  the  MPL,  keeping  its 
overall  porosity  constant  ~65%.  Two  parameters  are  varied,  the  ini¬ 
tial  pore  size  dpore  and  their  volume  fraction  r)p0re-  We  investigate 
the  following  five  cases: 

1 .  Case  I:  dPore  =  80  nm,  ^Pore  =  80% 

2.  Case  II:  dpore  =  160  nm,  t]p0Ye  =  40%. 

3.  Case  III:  dpore  =  160  nm,  r]p0re  =  20%. 

4.  Case  IV:  dpore  =  80  nm,  ^pore  =  20%. 

5.  Case  V:  dp ore  =  80  nm,  tjp ore  =  40%. 


for  the  concentration  c  in  the  pore  space.  This  solution  gives  us  the 
total  diffusion  flux  j  and  with  Fields  law 

(10) 


The  pore  size  distribution  of  these  models  is  given  in  Fig.  6. 
Increasing  the  initial  diameter  of  the  pores  (i.e.  Cases  II  and  III) 
results  in  a  more  densely  packed  structure.  This  is  observed  by  the 
formation  of  small  pores.  Decreasing  the  initial  solid  fraction  results 


j  —  ^bulk^C 
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Table  1 

Effect  of  structure  on  the  thermal  conductivity  and  diffusion  coefficient  of  the  MPL. 


Case 

Thermal 

conductivity 

Diffusion  coefficient 

(xlO-6  m2  s_1) 

Diffusivity 

(Q  =  ^)a 

I 

0.094 

1.92 

0.106 

II 

0.120 

1.78 

0.099 

III 

0.104 

1.65 

0.092 

IV 

0.092 

1.59 

0.088 

V 

0.095 

1.67 

0.093 

a  The  bulk  diffusion  coefficient  is  estimated  at  0°C  to  be  1.8  x  10-5  m2  s_1. 


in  shifting  the  pore  distribution  to  the  left;  in  other  words,  the  over¬ 
all  size  of  the  pores  decreases.  The  effect  of  these  two  parameters 
on  the  thermal  conductivity  and  diffusion  coefficient  are  calculated 
as  given  in  Table  1.  The  results  show  that  the  effect  of  the  struc¬ 
ture  on  the  overall  thermal  conductivity  and  diffusion  coefficient 
is  minimal.  For  Cases  II  and  III,  the  thermal  conductivity  is  slightly 
higher  than  the  base  case  (Case  I)  due  to  the  more  packed  agglomer¬ 
ates.  Similarly,  their  diffusion  coefficient  is  lower  due  to  the  higher 
contribution  of  the  Knudsen  diffusion  to  the  overall  diffusion  coef¬ 
ficient.  The  smaller  pore  sizes  are  also  the  main  contributor  to  the 
decrease  of  the  overall  diffusion  coefficient  of  Cases  IV  and  V. 

In  the  next  sections  of  this  paper,  the  MPL  under  investigation 
is  that  described  by  Case  I.  The  pore  distribution  of  this  structure 
and  its  diffusion  coefficient  compare  very  well  with  experimental 
data. 


4.2.  Pore  distribution  of  the  microporous  layer 

The  pore  distribution  of  the  reconstructed  MPL  is  given  in  Fig.  7 
and  compared  against  published  experimental  data.  As  it  can  be 
seen,  the  overall  trend  of  the  pore  distribution  is  similar  in  both 
cases.  The  radius  of  the  pores  of  the  reconstructed  MPL  lies  in  the 
range  0.01  <r<  0.085  pun  with  the  peak  occurring  at  a  radius  of 
0.04  p,m.  The  main  discrepancy  between  the  numerical  and  exper¬ 
imental  values  is  the  actual  height  of  the  curve.  This  is  explained 
by  the  fact  that  the  experimental  data  was  obtained  for  an  MPL 
applied  onto  carbon  paper  and  normalized  using  the  total  weight 
of  the  GDL/MPL  assembly,  which  is  much  higher  than  that  of  the 
MPL.  The  numerical  results  on  the  other  hand  are  obtained  for  the 
MPL  only. 


4.3.  Properties  of  the  microporous  layer 

For  the  purpose  of  simulating  the  transport  phenomena  in  PEM 
fuel  cells,  the  MPL  is  often  treated  as  a  stand-alone  domain  and  all 
the  conservation  equations  must  be  resolved  in  this  domain.  Using 
the  macro-homogeneous  assumption,  the  solid  and  void  regions  of 
this  layer  are  lumped  together  and  the  so-called  effective  properties 
are  used.  In  the  current  literature,  the  Bruggeman  approximation  is 
often  utilized  to  estimate  the  diffusion  coefficient  and  the  thermal 
conductivity  of  this  layer.  However,  as  illustrated  in  Figs.  8  and  9, 
this  approximation  does  not  accurately  predict  these  properties. 
The  Bruggeman  Approximation  overestimates  the  diffusion  coef¬ 
ficient  of  the  MPL  due  to  the  fact  that  the  Bruggeman  formula 
neglects  the  Knudsen  effect.  However,  the  Knudsen  effect  should 
not  be  neglected  when  the  pore  size  is  as  small  as  one  microm¬ 
eter  and  it  usually  becomes  dominant  when  the  pore  size  is  less 
than  1 00  nm,  which  is  the  case  for  the  MPL.  The  difference  between 
the  numerical  results  of  the  thermal  conductivity  and  those  pre¬ 
dicted  by  the  Bruggeman  Approximation  arises  from  the  difference 
in  structure.  The  Bruggeman  formula  was  derived  for  closely  packed 
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Fig.  7.  Comparison  of  the  pore  distribution  of  the  reconstructed  MPL  with  experi¬ 
mental  data  -  the  porosity  of  the  reconstructed  MPL  is  70%. 
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spherical  particles  whereas  in  the  case  of  the  MPL,  neighboring 
agglomerates  might  not  have  perfect  contact.  In  other  words,  the 
contact  resistance  between  the  spherical  agglomerates  is  different 
in  this  study  than  that  considered  by  Bruggeman  [30]. 

As  the  porosity  of  the  microporous  layer  directly  depends  on 
the  packing  of  the  agglomerates  and  the  MPL  is  usually  attached 
to  the  carbon  paper,  obtaining  its  physical  porosity  could  be  diffi¬ 
cult.  Experimentally,  two  porosity  values  have  been  reported  [20,9]. 
These  values  are  estimated  values  based  on  the  estimated  penetra¬ 
tion  and  thickness  of  the  layer  and  the  pore  volume  of  the  substrate. 
Gostick  et  al.  [20]  suggested  that  a  correlation  between  the  pore 
volume  and  the  thickness  of  the  layers  can  be  used  to  estimate  the 
porosity  of  the  MPL  as  follows: 

£mpl  =  (l-0)Gr-— (13) 
Cmpl 

where  6Mpl  is  the  porosity  of  the  MPL,  0  is  the  proportion  of  the  total 
pore  volume  in  the  portion  of  the  electrode  not  intruded  by  the  MPL 
material  and  is  found  using  curve  fits  of  the  capillary  pressure,  67 
is  the  total  porosity  of  the  electrode,  tT  is  the  total  thickness  of  the 
electrode  and  tMPL  is  the  thickness  of  the  MPL. 

Using  Eq.  (13),  Gostick  et  al.  [20]  estimated  the  porosity  of 
an  MPL  on  an  SGL  10BB  to  be  72%  and  Chan  et  al.  [9]  estimated 
the  porosity  of  an  MPL  applied  onto  SolviCore  carbon  paper  to 
be  64%.  The  variation  of  the  porosity  reported  by  the  two  studies 
could  be  associated  with  experimental  errors  as  well  as  differ¬ 
ences  in  structure.  As  mentioned  earlier  in  the  introduction  section, 
Chan  et  al.  [9],  measured  the  diffusion  coefficient  of  the  MPL 
to  be  1.54  ±0.21  x  10-6m2s-1,  which  translates  to  a  diffusivity, 
Q,  (Q.=  Deff/Dbulk)  of  0.073  ±0.01.  At  a  porosity  value  of  65%,  the 
numerical  estimation  of  the  diffusivity  of  an  oxygen-nitrogen  mix¬ 
ture  is  0.106.  The  difference  between  the  numerical  estimation 
of  this  study  and  the  experimentally  measured  diffusion  coeffi¬ 
cient  can  be  attributed  to  the  method  of  estimation  utilized  in  [9]. 
As  mentioned  earlier,  the  MPL  cannot  be  found  as  a  stand  alone 
layer;  hence,  the  diffusion  coefficient  measurements  are  done  using 
a  GDL/MPL  assembly.  In  [9],  the  authors  used  a  series  network 
resistance  to  estimate  the  diffusion  coefficient  of  the  MPL.  Due  to 
uncertainties  surrounding  the  penetration  of  the  MPL  into  this  GDL, 
this  method  does  not  take  into  account  the  resistance  to  diffusion 
due  to  this  penetration.  Hence,  the  diffusion  coefficient  of  the  MPL 
can  be  underestimated. 

As  mentioned  in  the  introduction  section  of  this  study,  the  ther¬ 
mal  conductivity  of  the  MPL  has  been  measured  experimentally 
and  discrepancy  between  these  measurements  is  apparent.  This 
discrepancy  can  be  the  result  of  various  reasons  with  the  most  sig¬ 
nificant  being  differences  in  the  actual  contact  area  between  the 
spherical  agglomerates.  In  the  MPL,  micro-cracks  could  be  present 
due  to  the  manufacturing  process.  These  cracks  can  result  in  a 
decrease  in  contact  area;  hence,  an  increase  in  thermal  resistance 
and  a  decrease  in  the  overall  thermal  conductivity.  Further,  exam¬ 
ining  the  thermal  conductivity  (of  a  SolviCore  MPL)  measured  by 
Burheim  et  al.  [13]  and  Unsworth  et  al.  [14]  another  reason  for 
this  difference  can  be  speculated.  This  difference  can  be  attributed 
to  the  amount  of  PTFE  added  in  the  production  of  the  MPL.  Exact 
information  in  regards  to  this  amount  is  not  available.  In  their  study 
Unsworth  et  al.  [14]  estimated  the  amount  of  PTFE  in  the  MPL  to  be 
approximately  30%. 

In  order  to  compare  the  numerical  results  of  this  study  to  the 
experimental  measurements  of  the  thermal  conductivity  of  the 
MPL,  one  needs  the  bulk  thermal  conductivity  of  the  carbon  mate¬ 
rial  used  for  the  manufacturing  of  the  MPL.  The  assumption  that 
the  carbon  of  the  microporous  layer  is  the  same  as  that  of  the  GDL 
could  be  questionable.  This  is  so  since  the  properties  of  the  carbon 
are  very  much  affected  by  the  heat  treatment  of  this  layer  dur¬ 
ing  manufacturing.  Further,  the  addition  of  PTFE  can  significantly 


reduce  the  bulk  thermal  conductivity  of  the  agglomerates.  The  ther¬ 
mal  conductivity  of  the  MPL  is  3  orders  of  magnitude  lower  than 
that  of  carbon  (see  [31,32]).  Comparing  our  numerical  result  in  this 
study  to  the  measurements  of  [13,14],  the  bulk  conductivity  would 
be  estimated  in  the  range  of  3-6WK-1  m. 


4.4.  Properties  of  the  GDL/MPL  assembly 

In  this  section,  the  effect  of  the  MPL  on  the  overall  properties 
of  the  electrode  is  investigated.  As  mentioned  earlier,  the  prop¬ 
erties  of  the  MPL  are  first  resolved  and  then  the  MPL  is  added  as  a 
homogenized  layer  in  the  electrode  and  the  overall  diffusion  coeffi¬ 
cient  and  thermal  conductivity  of  the  electrode  are  then  estimated. 
In  this  section,  we  investigate  the  effect  of  thickness,  porosity  and 
penetration  of  the  MPL  on  the  properties.  The  properties  used  for 
the  MPL  under  this  section  are  those  of  Case  I  described  earlier. 

The  penetration  of  the  MPL  into  the  GDL  is  achieved  by  filling  in 
the  pore  space  of  the  GDL  with  the  MPL  material.  In  this  case,  the 
effective  properties  calculated  for  the  MPL  using  its  model  are  used 
to  characterize  the  homogenized  MPL  material. 

The  penetration  and  thickness  of  the  MPL  are  highly  depen¬ 
dent  on  the  manufacturing  process  of  the  MPL  and  the  substrate 
on  which  it  is  applied.  From  experimental  observations,  these  two 
parameters  can  highly  vary.  In  Ref.  [9],  Chan  et  al.  measured  a  vary¬ 
ing  thickness  for  the  MPL  depending  on  the  substrate  on  which  it  is 
applied.  An  average  thickness  of  38.4  and  60.7  p,m  was  measured 
for  an  MPL  on  SIGRACET  25  DC  and  SolviCore  Type  B,  respectively. 
Further  in  Ref.  [20],  the  thickness  of  the  MPL  (applied  on  an  SGL 
carbon  paper)  has  been  estimated  at  about  118  pun  with  a  pene¬ 
tration  depth  of  63  [xm.  In  the  following  sections  of  this  paper,  our 
main  objective  is  to  numerically  understand  the  effects  of  these  two 
parameters. 


4.4 A.  Effect  of  thickness  of  the  microporous  layer 

As  mentioned  earlier,  the  thickness  of  the  microporous  layer  can 
vary  tremendously;  hence,  it  is  interesting  to  investigate  the  effect 
of  the  thickness  on  the  overall  diffusion  coefficient  and  thermal 
conductivity  of  the  electrode.  To  carry  out  this  investigation,  the 
MPL  porosity  is  kept  constant  at  65%  for  all  thicknesses.  The  effect 
of  the  MPL  thickness  on  the  diffusion  coefficient  and  the  thermal 
conductivity  is  illustrated  in  Fig.  10.  Due  to  the  small  diffusion  coef¬ 
ficient  of  the  MPL  in  comparison  to  the  GDL,  the  overall  diffusion 
coefficient  of  the  electrode  decreases  in  both  directions  (through- 
and  in-plane).  This  implies  that  the  small  pores  of  the  MPL  create 
higher  resistance  for  the  diffusion  and  increasing  the  thickness  of 
this  layer  will  result  in  hindering  the  overall  diffusion  process.  From 
Fig.  10,  the  effective  diffusion  coefficient  of  the  electrode  decreases 
with  the  increase  of  the  MPL  thickness.  This  implies  that  the  resis¬ 
tance  to  diffusion  due  to  the  addition  of  the  MPL  increases  with  its 
thickness.  However,  this  trend  is  not  valid  for  the  thermal  conduc¬ 
tivity.  As  it  can  be  seen,  the  through-plane  thermal  conductivity 
increases  with  the  increase  of  the  thickness  whereas  its  in-plane 
counterpart  decreases.  This  is  mainly  due  to  the  difference  between 
the  through-  and  in-  plane  thermal  conductivities  of  the  GDL.  The 
through-plane  thermal  conductivity  of  the  GDL  is  much  lower  than 
that  of  the  MPL;  hence,  the  addition  of  the  MPL  will  enhance  heat 
conduction  in  the  through-plane  direction.  In  the  in-plane  direc¬ 
tion,  the  thermal  conductivity  of  the  GDL  is  much  higher  than  the 
MPL;  hence,  the  addition  of  the  MPL  will  hinder  heat  conduction  in 
the  in-plane  direction. 

The  observed  trend  can  be  explained  from  a  resistance  net¬ 
work  perspective.  In  the  through-plane  direction,  a  series  network 
governs  heat  conduction  and  mass  diffusion  whereas  a  parallel  net¬ 
work  can  be  used  to  analyze  this  transfer  in  the  in-plane  direction. 
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Fig.  10.  Effect  of  the  thickness  of  the  MPL  on  the  overall  (a)  diffusion  coefficient;  (b)  Fig.  11.  Effect  of  the  porosity  of  the  MPL  on  the  overall  (a)  diffusion  coefficient;  (b) 
thermal  conductivity  of  the  electrode.  thermal  conductivity  of  the  electrode. 


The  resistance  through  each  layer  (GDL  and  MPL)  can  be  written 
as: 


(14) 


where  L  is  the  length  scale  through  which  the  transfer  occurs,  £ 
denotes  the  transport  property  and  Ac  is  the  cross-section  area. 

With  this  type  of  analysis,  the  through-  and  in-plane  properties 
of  the  electrode  can  be  obtained  as: 


£total  —  < 


_ ^total _ 

(f'MPL/ §MPL)  +  (f'GDL/ ^GDL) 
^MPL^MPL  +  §GDL^GDL 


L 


total 


Through-plane 


In-plane 


(15) 


where  the  superscripts,  total,  MPL  and  GDL  stand  for  the  electrode, 
microporous  layer  and  gas  diffusion  layer,  respectively. 

Fig.  10  also  shows  that  the  overall  thermal  conductivity  of  the 
electrode  is  higher  than  that  of  the  GDL,  whereas  the  diffusion  coef¬ 
ficient  is  lower.  This  is  the  case  since  the  thermal  conductivity  of 
the  MPL  is  higher  than  that  of  the  GDL  and  its  diffusion  coeffi¬ 
cient  is  much  lower.  Hence,  using  this  finding  and  Eq.  (15),  certain 
properties  of  the  electrode  under  compression  can  be  deduced.  As 


discussed  earlier  in  many  studies  [1  ],  the  MPL  is  very  mechanically 
stable.  Due  to  its  tightly  packed  agglomerates,  it  does  not  deform 
under  compression.  Hence,  under  compression  of  the  thermal  con¬ 
ductivity  of  the  electrode  (GDL  +  MPL)  would  be  lower  than  that  of 
the  bare  GDL.  This  phenomenon  is  also  easily  observed  from  the 
experimental  data  reported  by  Karimi  et  al.  [1 1  ].  Similarly,  the  dif¬ 
fusion  coefficient  of  the  electrode  can  become  higher  than  that  of 
the  bare  GDL  after  certain  compression  pressures. 


4.4.2.  Effect  of  porosity  of  the  microporous  layer 

The  porosity  is  a  measure  of  the  packing,  or  density,  of  the  micro¬ 
porous  layer.  It  is  interesting  to  understand  how  the  density  of  these 
agglomerates  would  affect  the  transport  coefficients  of  the  overall 
backing  layer  as  is  illustrated  in  Fig.  1 1 .  As  the  porosity  of  the  MPL 
increases,  the  diffusion  coefficient  of  the  backing  layer  increases 
whereas  the  overall  thermal  conductivity  decreases.  The  increase 
in  the  diffusion  coefficient  is  attributed  to  the  increased  pore  space 
in  the  MPL  available  for  diffusion.  At  the  same  time,  however,  due 
to  the  decrease  in  packing  density  of  the  MPL,  higher  thermal  resis¬ 
tances  are  introduced  and  the  thermal  conductivity  of  the  overall 
backing  layer  is  reduced. 
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Fig.  12.  Effect  of  the  penetration  depth  of  the  MPL  on  the  overall  (a)  diffusion 
coefficient;  (b)  thermal  conductivity  of  the  electrode. 

4.4.3.  Effect  of  MPL  penetration 

Due  to  the  small  size  of  the  carbon  agglomerates  (dia.  40  (nm)) 
of  the  MPL,  the  penetration  of  the  layer  into  the  carbon  paper  is 
inevitable.  The  average  pore  diameter  of  the  carbon  paper  GDL 
is  much  larger  than  the  agglomerate  diameter;  hence,  the  carbon 
agglomerates  will  fill  the  pores  of  the  GDL.  The  penetration  depth 
of  the  MPL  into  the  GDL  relies  heavily  on  the  manufacturing  pro¬ 
cess  of  the  electrode.  Due  to  this  dependence,  this  depth  varies 
tremendously  and  is  sometimes  difficult  to  accurately  characterize 
experimentally.  In  this  section,  we  study  the  effect  of  this  parameter 
from  a  numerical  perspective  to  gain  a  comprehensive  understand¬ 
ing  of  its  influence  on  the  diffusion  coefficient  and  the  thermal 
conductivity  of  the  porous  electrode.  The  dependence  of  these  two 
properties  is  found  to  be  linear  in  both  directions,  through-  and 
in-plane,  as  given  by  Fig.  12. 

5.  Conclusion 

Interest  in  the  microporous  layer  is  mainly  attributed  to  its 
ability  to  enhance  the  fuel  cell  performance.  However,  due  to  the 
complication  associated  with  the  separation  of  this  layer  from 
the  backing  layer,  experimental  measurements  of  its  transport 


properties  can  be  troublesome.  Using  stochastic  methods  to  recon¬ 
struct  the  3-D  pore  distribution,  we  were  able  to  investigate  the 
effect  of  the  pore  distribution  on  the  transport  properties  (mainly 
diffusion  coefficient  and  thermal  conductivity)  of  the  MPL.  The  esti¬ 
mated  thermal  conductivity  and  diffusion  coefficient  of  the  MPL 
were  used  to  carry  out  a  parametric  study  to  investigate  the  effect  of 
MPL  porosity,  thickness  and  penetration  depth  on  the  overall  prop¬ 
erties  of  an  MPL/GDL  assembly.  Further,  a  detailed  analysis  of  the 
differences  between  the  numerically  estimated  properties  of  the 
MPL  and  published  experimental  measurements  was  put  forward. 
Direct  comparison  between  the  numerical  findings  and  published 
experimental  data  is  quite  difficult.  The  numerically  estimated  dif¬ 
fusion  coefficient  was  found  to  be  higher  than  that  of  experimental 
measurements.  As  experimental  measurements  neglect  the  pene¬ 
tration  of  the  MPL  into  the  GDL,  they  tend  to  underestimate  the 
diffusion  coefficient.  The  diffusibility  of  an  MPL  with  a  65%  poros¬ 
ity  was  estimated  numerically  as  0.106.  In  the  case  of  the  thermal 
conductivity,  the  numerical  estimation  of  the  normalized  thermal 
conductivity  along  with  the  published  experimental  data  were  used 
to  estimate  the  bulk  thermal  conductivity  of  the  carbon  agglomer¬ 
ates  to  be  in  the  range  of  3-6WK-1  m. 
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